ag.state <- ag %>%
	group_by(state) %>%
	summarise(count.sum = length(state))

ag.state[49,"state"] <-"AK"
ag.state[50,"state"] <-"IL"
ag.state$state[ag.state$state == "NB"] <- "NE"

ag.state$count.sum[is.na(ag.state$count.sum)] <-0

(map <- plot_usmap(data = ag.state, values = "count.sum", color = "black") + 
		scale_fill_continuous(low = "white", high = "red", name = "Number of Sampled Agencies", label = scales::comma) + 
		theme(legend.position = c("right")))

ggsave(plot=map, "results/figure1.pdf", width=6,height=3)